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ABSTRACT 



We present the observational results of a survey designed to target and de- 
tect asteroids whose colors are similar to those of Vesta family members and 
thus may be considered as candidates for having a basaltic composition. Fifty 
basaltic candidates were selected with orbital elements that lie outside of the 
Vesta dynamical family. Optical and near-infrared spectra were used to assign 
a taxonomic type to 11 of the 50 candidates. Ten of these were spectroscop- 
ically confirmed as V-type asteroids, suggesting that most of the candidates 
are basaltic and can be used to constrain the distribution of basaltic material 
in the Main Belt. 

Using our catalog of V-type candidates and the success rate of the survey, we 
calculate unbiased size-frequency and semi-major axis distributions of V-type 
asteroids. These distributions, in addition to an estimate for the total mass 
of basaltic material, suggest that Vesta was the predominant contributor to 
the basaltic asteroid inventory of the Main Belt, however scattered planetes- 
imals from the inner Solar System (a < 2.0 AU) and other partially /fully 
differentiated bodies likely contributed to this inventory. In particular, we in- 
fer the presence of basaltic fragments in the vicinity of asteroid 15 Eunomia, 
which may be derived from a differentiated parent body in the middle Main 
Belt (2.5 < a < 2.8). We find no asteroidal evidence for a large number of 
previously undiscovered basaltic asteroids, which agrees with previous theo- 
ries suggesting that basaltic fragments from the ~ 100 differentiated parent 
bodies represented in meteorite collections have been "battered to bits" [Bur- 
bine, T.H., Meibom, A., Binzel, R.P., 1996. Mantle material in the Main Belt: 
Battered to bits? Met. k Planet. Sci. 31, 607]. 
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1 Introduction 



The surface properties of asteroids as determined by spectroscopic observa- 
tions are clues to the environmental conditions present in the early Solar Sys- 
tem. Asteroids with basaltic material on their surface are interpreted as trac- 
ers of bodies whose interiors reac hed the melting temperatur e of silicate rock 



and subsequently differentiated (IGaffey et al.l . Il993b , 



20021). Until recently 



Vesta and the dynamical 



observed surface basalt ( McCord et al. 



ly associated Vestoids were the only aster oids with 



1970 



Binzel and Xu 



19931 ). A small 



number of basaltic asteroids have now been f ound t hat may derive from other 



differentiated parent bodies ( 



Binzel et al. 



2007 



Roig et al. 



azzaro et al 



2000 



Hammergren et al 



2006 



20081 ). Other asteroids an d asteroid fa mihes 



may al so represent differentiated parent bodie s [e.g. Flora (IGaffevl. 



nomia (jReed 



at el. 



1997 



Nathues et al 



2005), Hebe flGaffev and Gilbert 



1984h. Eu- 



20051 ) ■ Eos (iMothe-Diniz and Carvanol. 



19981 ). and Thetis, Agnia and Merxia ( iSunshine et al. 



20041 )]. however it is not clear why so little basaltic material from these ob- 
jects remains while the Vestoids are one of the most prominent families in the 
Main Belt. Two possibilities are that the Vestoids were produced relatively 
recently and that non-Vestoid basaltic frag ments were heavily comminuted to 



sizes below our current observational 



1996 



Davis et al. 



2002 



Bottke et al. 



imits ( 



'. Davis et al. 



1989 



Burbine et al. 



2005bl ) or that some space weathering 



process has erased the diagno stic spectral signatures of surface basalts (e.g. 



Wetherill and Chapman 



19881) 



In addition to asteroidal evidence for the widespread occurrence of differ- 
entiation, oxygen isotope analyses of the iron meteorit es suggests that as 



many as ~ 60 differentiated parent bodies once existed (jChabot and Haack 
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2006 




Burbine et al. 


(Yang et al. 




2007 


)■ 



20071 ). Furthermore, the HED basahic achondrites are not iso- 



topica 



ly homogenous and ma y represent more than one differentiated parent 



body (lYamaguchi et al 



200l|). 



These populations suggest that basaltic material was common throughout 
the inner Solar System, however it is unclear whether there currently exists 
any basaltic vestige of the numerous differentiated parent bodies that once 
existed. To provide insight on these issues we have determined an unbiased 
distribution of V-type asteroids in the Main Belt. Other studies have calcu- 



lated unbiased dist r ibutions of asteroid typ es (e.g. 



Gradie et al. 



Gradie and Tedesco 



1982 



1989 



Bus and Binzel 



20021 ). but they have been restricted to 



the largest and brightest asteroids and did not consider the presence of V-type 
asteroids beyond 2.5 AU. Our survey provides constraints on the number of 
basaltic asteroids up to an absolute magnitude of 17.2, which is an increase 
of roughly 3 orders of absolute magnitude compared to these previous studies 
and therefore probes for basaltic fragments that may represent the upper size 
limit to a heavily comminuted population. 

Throughout this work we use specific terms related to the taxonomic classi- 
fication, surface mineralogy and dynamical properties of aste roids. V-type is 



a taxo nomic distinction based on spectral criteria defined by 



Bus and Binzel 



( 120021 ). Any object with photometric colors consistent with a V-type desig- 
nation is a V-type candidate. Basaltic asteroids are differentiated bodies or 
fragments of differentiated bodies with basalt on at least part of their sur- 
face. Basaltic asteroids are characterized based on interpreted mineralogy from 
spectroscopic data. An object with photometric colors indicative of a basaltic 



composition is a basaltic candidate. For this work we consider a V-type tax- 
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onomy to be directly attributable to a basaltic composition, therefore the 
terms V-type and basaltic are interchangeable. Vestoids are a dynamically- 
defined sub-population of V-types thought to have origin a ted fr om the Vesta 



19961 ) ■ Both spec- 



parent body at least 1 billion years ago (iMarzari et al.l . 
troscopic data and dynamical studies are required to classify an object as a 
Vestoid. Vestoids are typically interpreted as having a basaltic surface compo- 
sition, however some asteroids with basaltic material on their surface are not 
Vestoids, e.g. Magnya (j] 



Lazzaro et al. 



2OOOI ). An asteroid that is dynamically 



associated with Vesta but has not been spectroscopically observed is called a 
dynamical Vestoid. 

In the early Solar System the likely heat source for metal-silicate differentiation 



and the production of b asalt was the decay of s^ 



lort- 



such as ^^Al and '^"Fe 



601 



2000 



Grimm and McSweeru . 



Tachibana and Husd . 



1993 



i yed radioactive isotopes 



Goswami and Vanhala 



20031 ). The half-lives of these isotopes are 0.73 and 
1.5 Myr respectively, thus under this scenario a differentiated planetesimal or 
embryo must have accreted within a few half-lives (< 3 Myr after the onset of 
decay) to reach silicate melting temperatures (> 1000 K). An alternative heat 



source for material in the early Solar Syste m may have been e^ 



induc tive heating by an intense solar wind (IHerbert et al 



1991 



ectromagnetic 



Gaffev et al. 



2OO2I), however this mechanism is poorly constrained flMcSAveen et al.l. 120021) 



and m ay be incompatible with observations of T Tauri winds (IWood and Pellas 



199l|). 



The increase of accretion time-scale with heliocentric distance implies that 
planetesimals from the terrestrial planet region of the Solar System incor- 
porated a high abundance of live radionuclides and may have differentiated, 
whereas planetesimals that formed in the vicinity of the gas giants accreted 
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too slowly to h ave melted. It is unc ertain where the boundary between these 



two regions lay. 



Bottke et al. 



(120061) suggest that differentiation was primarily 
confined to the inner Solar System (a < 2.0 AU) so that the iron meteorites 
that reach the Earth today trace differentiated planetesimals that were scat- 
tered into the Main Belt from smaller heliocentric distances. It is unclear 
whether the presence of non-Vestoid basaltic aste roids and partially melted 



asteroids in the outer be 



(12005h : 



Reed at el. 



fll997h 



t [e.g. Eos and 



Nathues et al. 



vunomia, 



Mothe-Diniz and Carvano 



(120051 )] provide evidence for the in- 



situ different i ation of material in the Main Belt or can be explained by the 



Bottke et al 



(120061 ) scenario. 



Basaltic melts derived from a chondritic precursor will contain pyroxenes, 
olivine and plagioclase. Pyroxenes are spectroscopically characterized by deep 
and broad 1- and 2-fim absorption features. Figure 1 shows optical and near- 
infrared spectra of basaltic asteroids 38070 (1999 GG2) and 33881 (2000 
JK66), which are two of our basaltic cand idates (§2). The visible spectrum 



was obtained with ESI (Sheinis et al. 



with SpeX ( iRayner et al. 



2OO2I ) on Keck II and the near-infrared 



20031 ) on IRTF. Also plotted is a RELAB (Re- 
flectance Experiment Laboratory) spectrum of a basaltic eucrite meteorite. 
These objects clearly show the absorption features at 1 and 2 fim. The depth, 
center and shape of t hese features are dictated by temperature, grain size 



and calcium content (IHardersen et al. 



2OO4I ) and ar e attributable to Fe^ + 



ionic transitions within the pyroxene crystal structure (ISunshine et al. 



20041 ). 



Olivine and plagioclase also show Fe^"^ absorption features around 1 fim. 

Although the spectral features of basaltic material dominate in the near- 
infrared, optical data (0.4 — 0.95 fim) are easier to obtain for faint asteroids 
and can be diagnostic of such surface material. Thus, we primarily relied on 
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optical taxonomic classification as a way of selecting and characterizing our 
target asteroids. 



r 



Recent observational and theoretical work 



2006 



Bottke et al 



2006 



Binzel et al. 



2007 



jazzaro et al 



2000 



Nesvorny et al 



Hammergren et al. 



2008 



Roig et al. 



20081 ) along with the data presented here suggests that there exists non-Vestoid 
basaltic asteroids throughout the Main Belt that have yet to be discovered. 
We have placed new constraints on the overall population of these rare ob- 
jects by investigating relatively faint asteroids that are not typically included 
in spectroscopic surveys. One of the goals of all research related to the extent 
of thermal alteration in the early Solar System is to understand how many 
differentiated parent bodies are represented by the Main Belt asteroids. This 
work is a small step towards that understanding. 

The remaining sections of this paper are organized as follows: in §2 we outline 
the procedure of selecting target asteroids for our observational campaign. In 
§3 and §4 we discuss our observational procedure and data reduction tech- 
niques. In §5 we use our observational results to derive a size-frequency distri- 
bution of basaltic asteroids in the Main Belt. We use this result to constrain 
the semi-major axis distribution and the total mass of basaltic asteroids. In 
56 we discuss and summarize these results. 



2 Selecting Basaltic Asteroid Candidates 



We selected candidate basaltic asteroids from the photometr i c dat a in the 



third release of the Sloan Digital Sky 



Object Catalog (MOC, 



Ivezic et al 



Survey (SDSS, 



York et al. 



2000) Moving 



20021 ). Our targets were chosen based 
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upon their orbital elements and a statistically significant similarity to the 
ugriz photometric colors of the Vestoid population. The ugriz bands span a 
wavelength range from approximately 0.35 — 0.90 fim. These five filters provide 
a crude optical spectrum that can indicate the presence of basaltic spectral 
features (Fig. 1). 

The 204,305 moving objects in the MOC were filtered down to ~40,000 by 
removing multiple observations of single bodies, obje cts not matched to known 



aster oids in the ASTORB orbital element database (jJuric et al. 



2002 



Bowel] 



20071 ). and objects with any photometric error larger than one magnitude. If 
this constraint on the photometric errors was more rigid, then the number of 
final candidates would have decreased. This would have reduced the statistical 
significance of our result and made it more difficult to coordinate observing 
runs with observable targets. 

The SDSS magnitudes published in the MOC are typically denoted by u'g'r'i'z'. 



A transform (from the SDSS website, http://www.sdss.org/dr6/ ) is applied to 



these magnitudes so that solar correction in the AB system can be performed: 

u = u' - 0.04 
g = ^' + 0.060 (^'-r'-0.53) 
r = r' + 0.035 (r' - - 0.21) 
i = i' + OMl (r'-i'-0.21) 
z = z' - 0.030 {i' - z' - 0.09) + 0.02 (1) 

Four color combinations were calculated [u — g, g ~r, r — i, i — z) and solar 
corrected for each object in our filtered version of the MOC. The adopted 
solar colors were u — g = 1.43, g — r = 0.44, r — i = 0.11 and i — z = 0.03 



and were subtracted from each color permutation. I n addition, two 



component colors were calculated for each asteroid (iNesvorny et al. 



jrinci pal 



20051): 
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PCI = 0.396 {u-g) + 0.553 {g - r) + 0.567 {g - i) + 0.465 {g - z) 
PC2 = -0.819 iu~g) + 0.017 - r) + 0.090 - i) + 0.567 {g - z){2) 



These linear combinations define a set of axes that ma ximize the separation 



between S- and C-type asteroids (iNesvorny et al. 



20051). 



Using a Hierarchical Clustering Method (HCM, iZappala et al. 



1990) and a ve- 



locity cutoff of 70 m/s we identified 5575 dynamical Vestoids from a 2005 ver- 



sion of the AstDys database (http: / /hamilton. dm. unipi.it/cgi-bin/astdys/ astibo ) 



This cutoff does not include the 365 m/s escape velocity (IBinzel and Xu 



19931 ) 



from Vesta's potential well. Of these dynamical Vestoids, 1051 are included in 
the MOC. Certainly some fraction of these are non-V-type interlopers, how- 
ever Fig. 2 shows that the colors of these 1051 objects are distributed with 
a low dispersion which suggests that the number of interlopers is small and 
that this population can be used to represent the photometric colors of V-type 
asteroids. 

We fit histograms of these colors with gaussian profiles (Fig. 2) to define 
color ranges corresponding to percentile groups within the dynamical Vestoid 
population. These ranges were defined symmetrically relative to the fitted 
centroid of the profiles. For instance, the 50*'* percentile of dynamical Vestoids 
in the MOC have u — g colors between 0.17 and 0.35 with a most probable 
value of 0.26. We used the 50*'' percentile for each of the 6 colors to statistically 
define the SDSS photometric properties of candidate basaltic asteroids. Table 
1 shows the adopted range for each color. Basaltic candidates were selected as 
objects whose colors overlap all 6 of these ranges. 



Dynamical constraints were then applied to objects that met the photomet- 
ric criteria. Asteroids were excluded from the candidate list if their orbital 
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properties overlapped the region of phase space occupied by the dynamical 
Vestoids. This dynamical constraint was applied to the osculating values of 
semi-major axis, eccentricity and inclination. Proper orbital elements would 
have been preferred but these are available for only a subset of MOC objects. 
The extent of phase space occupied by the 1051 dynamical Vestoids in the 
MOC is 2.22 - 2.48 AU in semi-major axis, 0.029 - 0.169 in eccentricities, 
and 4.82 - 8.26° in inclination. These photometric and dynamical constraints 
produced 50 targets for spectroscopic follow-up (Table 2, Fig. 3). 

Table 2 indicates the relative ejection velocity, Av, required for each of the 
candidates to have originated from the surface of Vesta. Simplified forms of 
Gauss's equations can be used to express an asteroid's change of velocity in 
three components (tangential, radial and perpend icular to the orbital p lane) 



Zappala et al. 



19961). The 



relative to another body (Equations 4, 5 and 6 in 
net relative ejection velocity, At>, is calculated by adopting these authors' 
values for the true anomaly and argument of perihelion of the Vesta family 
at the time of its collisional formation (see their Fig. 12) and then adding 
the three components in quadrature. The 365 m/s escape velocity from Vesta 



(iBinzel and Xu 



19931 ) is not taken into account for these calculations. With 
the exception of two objects [28517 (2000 DD7) and 156914 (2003 FB31)] all 
of the ejection velocities are greate r than the ~ 0. 6 km/s that is expected for 



multi-kilometer Vestoid fragments (lAphaug 



19971). 



It should be noted that calculating Av for a pair of asteroids can be a useful 
metric for quantifying dynamical proximity, however dynamical effects such 
as Yarkovsky drift and resonance scattering can produce values of Av for 
members of the same family that are much large r than the sub- km/s veloci- 



ties expected from HCM family definitions (e.g. 



Zappala et al. 



1995, 



19961) 
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and from SPH simul ations of fam ily forming impacts ( e.g. lAphaugi . 11997 



Nesvorny et al. 



20061 ). For insta nce, 



with Av > 2 km / s are likely and 



Nesvorny et 



Roig et al. 



al. 



(120081 ) show that Vestoids 



(j2008|) suggest that asteroid 40521 



(1999 RL95) is a fragment from Vesta's surface even though its Av = 2.3 km/ s. 
Part of this problem is rooted in the HCM technique: only objects that "stand 
out" from the background population of asteroids can be identified as mem- 
bers of a collisional family. In reality the dynamical boundaries of collisional 
families are not so clearly defined and certainly extend into the background 
population and even other families. 

For these reasons it is likely that our HCM-defined dynamical boundaries do 
not encompass the full extent of the Vestoid family, a presumption that is 
supported by the presence of basaltic candidates near these boundaries (Fig. 
3). However, the inclusion of these objects as candidates does not dramatically 
affect our analysis or conclusions. Using a larger cutoff velocity (> 70 m/s) 
to expand the dynamical boundaries would connect all of the asteroids in the 
inner Main Belt (a < 2.5 AU) and would therefore neglect the possibility of 
selecting basaltic asteroids from other Vesta-like differentiated parent bodies. 
The use of a "small" cutoff velocity (i.e. one that does not encompass the entire 
Vestoid population) is conservative in the sense that it makes the amount 
of supposedly non- Vestoid basaltic asteroids larger and thus contributes to 
raising the upper limit for the total mass of basaltic asteroids (§5.2). 



Other studies have employed si milar color-selection techniques for identifying 



V-type candidates in the MOC (IRoig and Gil-Hutton 



2006 



Binzel et al. 



20071). 



Roig and Gil-HuttonI (120061 ) used two principal com- 



2006 : 



Hammergren et al. 



ponents and i—z colors to place photometric constraints on the MOC, generat- 
ing over 250 non- Vestoid basaltic candidates. Sixteen of our 50 candidates are 
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included in the 



Roig and Gil-HuttonI (120061 ) list (Table 2). The lack of o verlap 



Roig and Gil-Huttonl (120061 ) used 



is due to the different selection techniques 
only three colors to select candidates, but applied more stringent criteria on 
the acceptable photometric error bars. Our techniq ue used 6 colors but applied 



more relaxed constraints on the photometric errors. IHammergren et al. 



(l2006h 



selected V-type candidates based on i — z co 
1 fim basaltic absorption feature (Fig. 1) and 



ors which are diag nostic of the 



Binzel et al. 



(120071 ) developed a 



system to translate the photometric measurements of MOC o 



j.iects into taxo- 



nomic classifications based on additional spectroscopic data (IBus and Binzel 



20021 ) . The last three columns in Table 2 indicate whether our candidates were 



also identified as V-type candidates by these other groups. 



It is unclear whether one of these methods i s more effective at taxonomi c pre- 



diction than the others. Our method a ndthelRoig and Gil-Hutton 



nique produce similar results: 11 of 12 



Roig and Gil-Huttonl (120061 ) candidates 



(I2OO6I ) tech- 



and 1 out of 11 of our t argets are spectroscop i cally 



work. IroIs et al 



(l2008h 



Roig and Gil-Huttonl (l2006f ): 



c onfirmed V-typ e s (thi s 



Lazzaro et al. 



(l2004h : 



Bus and Binzell (120021)'). The significant ove rlap between our candidates and 



those selected by 



Hammergren et a 



are roughly equivalent. 



Binzel et al. 



(120061) suggest that these two methods 



(120071 ) and Masi et al. (2008, submitted 



to A&A) show that their technique is ~ 80% effective at predicting taxonomic 
type. 



14 



3 Spectroscopic Observations 



We have obtained spectra of 9 basaltic candidates|_J Table 3 contains a sum- 
mary of the observational circumstances for each target asteroid. In addition, 
the candidates 40521 (1999 RL95) and 2823 (van der Laan) have been observed 



as part of comp 



2001 



Roig et al. 



ementary projects performed by other groups (IBinzel et al. 



20081 ) 



We u sed the Supernova Integral Field Spectrograph (SNIPS) fiLantz et al. 



2004J ) on the University of Hawaii 2.2 meter telescope to obtain optical spectra 
(0.4 — 95 fim) of asteroids brighter than V ~ 18.5 magnitude. The SNIPS 6" x 
6" field of view is divided into 0.4" x 0.4" spaxels (spectral pixels), which are 
each generated by one lens in a 225 element multi-lenslet array. The light from 
each spaxel is fiber-fed to a dichroic mirror that splits the light into red and 
blue channels. The two channels pass through separate grisms and are imaged 
on two 2k X 4k blue- and red-optimized CCDs. The approximate spectral 
resolution (R) of this instrument is 1000 at 0.43 fim and 1300 at 0.76 fim. Bias 
frames, dome fiats and internal arc lamp spectra are automatically obtained 
as part of the standard SNIPS observing routine. 



The Echellette Spectrograph and Imager (ESI) on Keck II (ISheinis et al. 



2OO2I ) was used to obtain optical spectra (0.38 — 0.90 fim) of objects that 



were too faint to observe with SNIPS. ESI was operated in the low resolution 



^ We have also observed a number of asteroids that fall just outside of the 50*'' 
percentile color rang e s (no t listed in Table 2) and asteroids that are part of the 



Roig and Gil-HuttonI (j2006l ) candidate list. The taxonomic results of these observa- 



tions can be found at www.ifa.hawaii.edu/~nmosko/Home/astobs.html 
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{R ~ 1000) prism mode and with a wide slit of 6". We found that this width 
resulted in greater consistency between multiple spectra of the same asteroid. 
To minimize the effect of atmospheric dispersion the orientation of the slit 
was set close to the parallactic angle at the midpoint of the observation. Bias 
frames and dome flat field images were taken at the beginning and end of each 
night. Wavelength calibration was achieved by taking spectra of internal Hg, 
Cu, and Xe arc lamps at the beginning and end of each night. The dispersion 
for ESI was found to be very stable so that one set of arc line images could 
have been used to calibrate all of our data. 



Near -infrared (NIR) observations were performed with SpeX (IRayner et al. 



20031 ) on NASA's Infrared Telescope Facility (IRTF). The telescope was op- 
erated in a standard ABBA nod pattern and SpeX was configured in its low 
resolution (R=250) prism mode with a 0.8" slit for wavelength coverage from 
0.8 - 2.5 /im. Internal fiat fields and arc line spectra were obtain ed immediately 



after each target asteroid. The IDL-based SpeXtool package ( jCushing et al 



20041 ) was used for data processing. 



All targets were observed at an air mass of less than 1.3 and as near as possible 
to their meridian crossings. Bright field stars were used to guide the telescopes 
at non-sidereal rates so that the target asteroid remained centered on the slit 
for ESI and on the lenslet array for SNIPS. For the optical instruments we 
employed a maximum exposure per image of 900 seconds. For SpeX, exposure 
times were limited to 120 seconds. The number of exposures was selected in 
an attempt to obtain a consistent S/N 100 for all asteroids. 

A solar analog star was observed to perform calibrations for each asteroid spec- 
tra (i.e. solar correction, compensating for instrument response and removal 
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of atmospheric absorption features, see §4). Analogs were observed closely 
in time, altitude and azimuth to the asteroid. Each one-dimensional asteroid 
spectrum was divided by that of an analog during the reduction process. For 
SNIPS and SpeX data a single solar analog was selected for each asteroid. We 
found that this prescription produced the most consistent results, which may 
be related to the use of brighter rat her than well es t ablish ed solar analogs. 
For ESI we followed the procedure of IBus and Binzell (120021 ) and divided the 



asteroid spectra by the average of all analogs for a given night. The specific 
analog that was used for the calibration of each asteroid is given in Table 3. 



4 Data Reduction 

Although the reduction of asteroid spectral data varies from one instrument 
to the next, a number of key steps are common to all optical reductions. 
We outline non-standard steps here because we are not aware of an earlier 
comprehensive description of thes e processes. For furth er details on spectral 



processing of our optical data see (IWillman et al 



20081) 



Raw two-dimensional spectral images have a number of superimposed com- 
ponents that must be removed to isolate one- dimensional asteroid reflectance 
spectra. Pixel-to-pixel variations on a CCD must be corrected by dividing data 
images with a normalized flat fleld image of a uniformly illuminated source. 
The normalization of this flat fleld must take into account the sensitivity vari- 
ations on the CCD as function of position and as a function of wavelength. 
Some CCDs display optical fringing at red wavelengths (the ESI CCD shows 
fringing effects at the few percent level around 0.9 /im), which can be cor- 
rected by dividing the data images with a master fringe frame. The fringe 
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pattern that affects data images is produced by parallel light rays originat- 
ing at infinity, thus it can not be reproduced with dome flats. A normalized 
exposure of the twilight sky can be used to reproduce fringe patterns. The 
extraction of the ID spectrum from the 2D image is performed by deflning 
and subtracting an average background from the asteroid data. This effec- 
tively removes any atmospheric emission features. To correct for atmospheric 
absorption the ID asteroid spectrum is divided by a solar analog spectrum. 
Division by the analog also removes the superimposed solar contribution in 
the measured reflectance. 

For SNIFS, one- dimensional spectra from the red and blue channels are pro- 
duced within a few minutes of data acquisition. These spectra are bias level- 
calibrated, flat fleld-corrected, cosmic ray-cleaned and have an assigned wave- 
length solution. We have developed a series of custom IDL scripts to perform 
the remainder of the data reduction. These scripts normalize the spectra at 
0.55 /im, combine multiple exposures of a single object, join the red and blue 
channel data, divide the asteroid by a solar analog and flnally re-bin the data 
to increase S/N. Around 0.5 /xm, the sensitivity of the SNIFS red and blue 
channels drops, thus causing signiflcant degradation of signal-to-noise. For 
cosmetic purposes the data in this region was removed. 

Spectral images from ESI have a rotated dispersion axis relative to the x-axis 
of the CCD and signiflcant spectral line curvature along the spatial direction. 
We corrected for these two effects by rotating the 2D images by 6.4° clock-wise 
and then straightening the curvature by applying a transform map produced 
from arc lin e images. The identify, reidentify, fitcoords and transform IRAF 



(Tody 



19861 ) routines were used to perform this straightening. Once all of the 



spectral images were straightened the reduction proceeded following standard 
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protocols. Dispersion solutions were determined based upon 28 different arc 
lines distributed throughout the wavelength range of the instrument. After 
dividing the asteroid spectrum by the solar analog we further divided it by a 
scaled mo (del of the atmospheri c absorption features that are present beyond 
0.85 nm jwillman et al.i . i2008l ) This procedure is a fine tuning to the final 
asteroid spectrum and effectively eliminates any residual absorption. 

Reduced optical spectra of the observed candidates are presented in Fig. 4. 
On the night of 18 June 2007, the blue channel spectrograph on SNIPS was 
inoperable, thus only the red channel data are presented for asteroid 60669 
(2000 GE4). A spectrum of Vestoid 6093 Makoto is included in Fig. 4 for 
reference; this asteroid is not one of our basaltic candidates. Only NIR obser- 
vations were obtained for candidate 33881 (2000 JK66). The SpeX spectrum 
of this target is shown in Fig. 1 and clearly shows the deep 1- and 2-micron 
absorption bands that are indicative of a basaltic surface composition. 

Asteroid 10537 (1991 RY16) was observed with both SNIFS and ESI (Fig. 5) 
and thus allows a comparison between spectra from the two instruments. ESI 
appears to have better red sensitivity, however we have not tried to fine tune 
the SNIFS data with an atmospheric absorption model. The slight discrep- 
ancy between these two spectra at the red-most wavelengths is likely due to 
imperfect correction of atmospheric features and the difference in solar analog 
selection. Other than these subtle differences there is very good consistency 
between the SNIFS and ESI data. 

A number of criteria were used to discriminate a V-type asteroid from the 



optical data (IBus and Binzel 



2002). These included: (1) a moderate .4 



0.7 nm slope that can include a shoulder at ~ 0.6 fim (Vilas et al 



20001), (2 
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a peak reflectance of ~ 1.2 at 0.75 /im wlien normalized at the canonical 0.55 
fim and (3) the blue edge of the 1 iim absorption feature reaching a depth of 
0.8 or greater at 0.9 /im. 



Asteroid types A, Q and R have similar spectral features to V-types and could 
become false positives in the taxonomic classification of our targets. However, 
errors produced by incorrectly classifying an A, Q or R-type asteroid have 
little bearing on our conclusions since it is l ikely that A- and R- t ypes are also 



deriv ed from differentiated parent bodies (IMarchi et al. 



2004 



Lucey et al. 



19981 ). Furthermore, since there are no Main Belt Q-types amongst the more 
than 2000 asteroids with assigned taxonomies it is improbable that any of our 
fifty Main Belt candidates are Q-type. 



In some cases a taxonomic distinction can be made (iBus and Binzell . |2002| ). A- 
types tend to have steeper slopes blue- ward of 0.75 fim and shallower 1.0 fim 
absorption features. At optical wavelengths Q-types are typically distinguished 
from the V-class by a smaller peak refiectance (~ 1.1) and a shallower 1.0 fim 
absorption feature. R-type asteroids have very similar spectral features to 
those of the V-class, however these asteroids are exceedingly rare. It would be 
unlikely that we discover one of these as part of our observational program. A 
distinction between R- and V-types can be made based on a shallower 1 iim 
absorption feature for R-types. 

S-complex asteroids are the most common contaminant in our list of candi- 
dates. This is attributable to a few factors. By number, they are the second 
most abundant asteroid type in the Main Belt. In addition, the difference 
between S and V-class asteroids can be obscured by the relatively low signal- 
to-noise of the MOC measurements and the statistical method by which we 



20 



have selected these targets. Fortunately S-types are easy to spectroscopically 
distinguish from V-types. For an S-type, the typical slope blue- ward of 0.75 /im 
is less steep and the 1 iim absorption feature is shallower (~ 0.9). 

The spectral classifications of our observed basaltic candidates are listed in 
Table 3. The truncated spectrum of 60669 (2000 GE4) was sufficient to de- 
termine its taxonomic class. The spectrum of 24941 (1997 JM14) is a limiting 
case for a V-type asteroid, if the relative reflectance at 0.9 /im were any larger 
then it would be classified as an S-type. The one candidate that turned out 
to be an S-type was 46690 (1997 AN23). Even with the low S/N of its spec- 
trum it is clear that the 1 /im absorption feature is not as deep as that of a 
typical V-type. As discussed in §6 there is some uncertainty about whether 
10537 (1991 RY16) is a basaltic asteroid, but we consider the observation of 
this target a successful confirmation of a V-type because of the depth of its 1 
/im band. 

5 Size and Orbital Distributions of Basaltic Asteroids 

Using the results and construct of our survey we develop a simple method- 
ology to determine unbiased number distributions of basaltic asteroids that 
are not part of the dynamical Vestoid population. We calculate a differential 
size-frequency distribution, N{H) where H is absolute magnitude (a proxy 
for size), and a semi- major axis distribution, N{a; H). We use N{H) to place 
estimates on the mass of basaltic material represented by our survey. These 
calculations are performed by correcting our candidate distributions for bi- 
ases inherent to the selection process and the SDSS catalog. The underlying 
assumption for these calculations is that the 11 basaltic candidates that have 



21 



been observed are taxonomically representative of the other 39. 



5.1 Size- Frequency Distribution 

We begin by equating the expected number of basaltic asteroids in our survey 
to the number represented by the observations: 



where N{H) is the unbiased, differential number distribution as a function 
of absolute magnitude for basaltic asteroids outside of the dynamical Vestoid 
population. This is the quantity for which we want to solve. Pg is the prob- 
abihty that an asteroid with basaltic surface material would be selected as a 
candidate based on its photometric properties. C{H) accounts for the com- 
pleteness and efficiency of the SDSS MOC. Completeness and efficiency are 
two separate properties. Completeness is the fraction of all asteroids detected 
by a given survey whereas efficiency is the fraction of observationally acces- 
sible asteroids that were detected. The MOC is not 100% complete for any 
value of H because it did not have full sky coverage. However, for objects 
brighter than H — 14.5 it is 100% efficient, because it detected all of the 
objects that were accessible to the survey up to that limiting magnitude. F 
is the false-positive rate of our survey. Since ten of eleven basaltic candidates 
have been spectroscopically classified as V-type asteroids (Table 3) the false 
positive rate of basaltic detection, F = 0.09. As stated in §2 this value may 
be elevated due to uncertainties in defining the boundaries of the dynamical 
Vestoid population. For instance, the inner Main Belt (a < 2.5 AU) candidates 
may be characterized by one value of F that reflects the scattering of Vestoids 




(3) 
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to large values of Av while another value characterizes the candidates beyond 
2.5 AU. Without additional observations of middle and outer belt candidates 
we adopt a single value of F for the entire Main Belt and therefore produce 
an upper limit to the basaltic material represented by the survey. The final 
parameter in Eq. 3, n{H), is the differential number distribution of our 50 
basaltic candidates. 

Ps is calculated as the fraction of the 1051 dynamical Vestoids in the MOC 
that would be selected as candidates based on their photometric colors (this 
inherently assumes that these dynamical Vestoids are composed of basaltic 
material). We find that Pg = 0.085, a low value that is probably related to our 
stringent constraints across 6 different colors and the relatively large errors on 
the SDSS photometry. 

C{H) is derived by comparing the SDSS MOC to the ASTORB database of all 
known asteroids. The ratio of the differential absolute magnitude distributions 
of these two databases gives the product of the efficiency and completeness 
of the MOC relative to ASTORB. However ASTORB is not 100% complete 
past H ~ 14.8. We extrapolate the number of asteroids for H > 14.8 with a 



2OO2I). 



function of the form k10°'^ where k and a are constants (jjedicke et al. 
Fitting the ASTORB distribution between H = 12.9 and H = 14.8 gives the 
unbiased number of asteroids as a function of absolute magnitude: 

NaiiiH) = (5.9 X 10-^) ■ 10°-^^. (4) 



We limit the extrapolation of this function to a narrow range beyond H = 14.8. 
The single slope approximation will slightly overestimate the population of 
basaltic asteroids because the slope of the distribution is known to decrease 
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at larger H (iJedicke et al.l . 120021 ) . 



Our 50 basaltic candidates span a range of absolute magnitudes from 12.9 
to 17.2. We focus all of our calculations on the number of basaltic asteroids 
within this range. Figure 6 plots the number distributions of the MOC and 
ASTORB between H = 12.9 and 17.2. The fit to the ASTORB data (Eq. 4) 
is indicated by the dashed line. The ratio of the MOC distribution to Nau{H) 
yields C{H). We find that the combined efficiency and completeness of the 
MOC is approximately 22% out to H = 14.5 and then drops off sharply to 
less than 1% at H = 16.5. 

We solve Eq. 3 for N{H) and place upper and lower 90% confidence levels on 
the calculated values (Fig. 7). The upper limits are equal to the means of the 
Poisson distributions for which there is a 90% chance of observing n{H) + 1 
candidates within a given bin, the lower limits are equal to the means for 
which there is a 90% chance of n{H) — 1 candidates. Only upper limits can 
be calculated for the if-bins with no candidates. The apparent fiattening of 
N{H) between if = 13 and 15 is likely a consequence of low number statistics. 

A single function of the form k10°^ was fit to the N{H) distribution by chi- 
squared minimization (Fig. 7). Extrapolation of this fit suggests that there 
are no basaltic asteroids with H < 11 outside of the dynamical Vestoid pop- 
ulation; a result that is confirmed (with the exception of Magnya) by spec- 
troscopic surveys which are close to 100 % complete at if = 10 and drop to a 



Bus and Binzell . 



2002 : 



Lazzaro et al. 



comp leteness of ~ 30% at ii = 11 (e.g 
20041 ). Combining our calculated N{H) distribution with the findings of these 
spectroscopic surveys suggests that the overall distribution of basaltic aster- 
oids outside of the dynamical Vestoid population is heavily weighted towards 
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small bodies < 5 km in size {H > 13). The implications of this result are 
discussed in §6. 

N{H) is used to calculate a cumulative number distribution of basaltic as- 
teroids (Fig. 8). To place an upper limit on the amount of basaltic material 
represented by our survey we set the cumulative number of asteroids within a 
given absolute magnitude range equal to the 90% Poissonian upper limit. In 
§5.3 the unbiased cumulative number of basaltic asteroids is used to normalize 
the semi-major axis distributions. 



5.2 Masses of Basaltic Material 



We estimate an upper limit to the mass of basaltic material represented in 
our survey by using the N{H) re sult. The diameter o f an asteroid, D, with an 



assumed geometric albedo, p, is (IBowell et al. 
1347 km 



19891): 



D 



-10 



-0.2H 



(5) 



Assuming constant density, p, and a spherical body we write Eq. 5 in terms 
of mass: 



M{H) = (1.28 X 10'* kg)^10 



p 



(6) 



Thus the total mass represented by our survey is given by: 

17.2 17.2 



: j M{H) N{H) dH = (1.5 X 10^2 kg) J N{H) iQ-o ^^ dH, (7) 

12.9 12.9 



where the integral is taken over the range of H represented by our candidates, p 
is set to 3000 kg/m^ and p is assumed to be 0.4, a value that is approximately 
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Vesta's albedo. We solve this integral by integrating over the N{H) upper 
limits and obtain Mgur = 5.9 x 10^^ kg. The use of the N{H) upper limits 
produces a liberal approximation to the total mass represented by our survey. 
We emphasize that this mass represents the amount of basaltic material within 
a limited ra nge of absolute magni t ude. Based on the r esults of spectroscopic 



surveys (e.g. 



Bus and Binzel 



2002 



Lazzaro et al 



20041 ) we expect there to be 



few undiscovered basaltic asteroids with H < 12.9, thus the lower H-limit to 
our calculation is not expected to significantly impact the mass result. Beyond 
the upper H-limit there may be enough basaltic material to increase Mg^r by 
a factor of a few depending on the number distribution at large H. 

To place M^^j, in perspective, we conservatively estimate the mass of basaltic 
material excavated from the surface of Vesta (based on the size of its largest 
crater) and place an upper limit on the mass of basaltic material represented 
by the dynamical Vestoids (My). As we show in §6 this upper limit to the 
total mass of basaltic material in the Main Belt [Mg^r + My) is still much 
less than that expected from the volume of Vesta's largest crater. 



The volume of Vesta's largest crater is estimate d to be 10^ km^ (IThomas et a. 



19971 ) and Vesta's mean density is 3.3 g/cm^ (jViateau and Rapaport 



This volume of material was excavated anywhere from 1-3.5 Gyr ago (IBottke et al. 



200lh. 



2005a 



Marzari et al. 



19961 ) and since then has been depleted by collisional 
grinding and the scattering of fragments out of the asteroid belt. Such deple- 
tion involves a complicated interplay of processes whose efficacy is dependent 
on the size distribution of ejected fragments. Det ailed collisional a nd dy nami- 



cal models are necessary to address this process. 



Nesvorny et al. 



(120081 ) show 



that ~ 10% of Vestoids would dynamic ally escape th e asteroid belt over a time 



scale of 2 Gyr. Based on the results of .Bottke et al 



(l2005bl ) we estimate that 
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50% is a reasonable estimate for the coUisional depletion of the observable 
Vestoid family (sizes > 1 km) over a lifetime of 3.5 Gyr. Thus we expect the 
combined collisional and dynamical depletion of the excavated crater mate- 
rial to be approximately 60%. Accounting for this depletion, we calculate the 
remaining mass of excavated basaltic material, M^.^ = 1.3 x 10^^ kg. 



It is unlikely that a significant fraction of Mi,^ is non-basaltic. The excavation 
depth, based o n the size of Vesta's largest crater, was approximately 13 km 



( jThomas et al. 



19971 ). which is less th an o r 



the thickness of Vesta's basaltic crust (IKeil 



comparable to calculations for 



size of t he largest basalt i c Vest oids, 1929 KoUaa (iKelley et al 



Peking (iBus and Binzel 



2002) 



20021) and is cornpara ble to the 



20031 ) and 2045 



We calculate the net mass of basaltic material in our list of 5575 dynamical 
Vestoids by de-biasing this population and removing non-V-type interlopers. 
Table 4 lists interlopers in the Vestoid dynamical space that are spectroscopi- 
cally confirmed to be non-V-type and have H < 13. Objects fainter than this 
typically do not have an assigned taxonomy, thus we assume that the fraction 
of interlopers for if > 13 is ~ 20%, the same as the percentage of interlopers 
for H < 13. The absolute magnitude distribution for dynamical Vestoids with 
non-V-type interlopers removed is shown in Fig. 9. 



Our list of dynamical Vestoids is derived from the AstDys database which is 
complete to an absolute magnitude of H = 14.8. We extrapolate the number 
of basaltic dynamical Vestoids for H > 14.8 using the same methodology as 
in §5.1. The fit to the dynamical Vestoid distribution from 13.4 < H < 14.8 
is shown as the dashed line in Fig. 9. The total basaltic mass in the dynami- 
cal Vestoid population. My, is calculated with Eq. 7 by replacing N{H) with 
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the differential number distribution of dynamical Vestoids, Nv{H), which is a 
combination of the interloper-corrected number of objects for Ny{H < 14.8) 
and the fitted number of objects for Ny{H > 14.8). The limits of integra- 
tion are changed to 11.9 and 18.0 so that all known dynamical Vestoids are 
included. The result is My = 4.8 x 10^*^ kg. 



A few issues may have slightly affected the calculation of My. First, a single- 
slope fit to the dynamical Vestoid distribution (dashed line in Fig. 9) likely 
produced an overestimate because unb iased absolute magn itude distributions 



are expected to turnover at large H fiJedicke et al. 



2OO2I ). Second, the use 



of a more recent version of the Ast Dys database wou . 



Vestoid family with more members (INesvorny et al. 



d pr o duce a dynamica l 



2008 



Roig et al. 



20081) 



and thus could slightly increase the completeness limit of this population. 
However, the additional dynamical Vestoids in an updated database would 
have H > 14.8 and therefore it is unlikely that they would greatly influence 
the mass estimate which is based on a reasonable extrapolation for the number 
of objects at values of if > 14.8. Lastly, My may be slightly overestimated 
due to unidentified interlopers in the dynamical Vestoid orbital space. All of 
these issues are expected to produce minor effects and as a whole are expected 
to increase the upper limit on My. We discuss the calculated masses, M^ur, 
Mex and My, in §6. 



5. 3 Semi-major Axis Distribution 

To calculate the semi-major axis distribution of basaltic asteroids we apply the 
same methodology as that used to derive the absolute magnitude distribution, 
i.e. we simply rewrite Eq. 3 in terms of the independent variable a, the semi- 
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major axis and H: 



N{a;H) ^ A{H) 



(1 - F) n{a) 
Ps C(a, H) ' 



(8) 



where A[H) is a normalization factor determined from the N{H) cumulative 
number distribution. A{H) is set so that the integrated number of basaltic as- 
teroids across all semi-major axes is in agreement with the absolute magnitude 
result. 

The completeness and efficiency of the MOC as a function of semi-major 
axis, C(a, H), is also dependent on absolute magnitude. For all of the basaltic 
asteroids represented by our survey: 



where nsoss is the observed semi-major axis distribution of asteroids in the 
SDSS MOC with H less than the specified limit and Nast is the unbiased semi- 
major axis distribution from the ASTORB catalog. However, the latter is only 
unbiased for H < 14.8. We calculate N^sTiO'] H < 17.2) by assuming it is a 
scaled version of the A^astI^; H < 14.8) distribution. The scaling is performed 
so that the integrated number of asteroids represented by NAsria; H < 17.2) 
is equal to the cumulative number predicted by Naii (Eq. 4). An identical 
procedure is followed to calculate C(a, H) for other absolute magnitude hmits. 

Figure 10 shows the unbiased semi-major axis distributions, N{a] H), for three 
different absolute magnitude limits. 90% Poissonian confidence levels are plot- 
ted as the upper and lower limits. Some of the semi-major axis bins have zero- 
value due to the scarcity of candidates with H < 15.0 and a > 2.5 AU. The 

apparent overlap of these non-Vestoid basaltic asteroids with the region occu- 



C(a; H < 17.2) 



nsDSsja; H < 17.2) 
NAsria; H < 17.2) ' 



(9) 
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pied by the dynamical Vestoids is a consequence of projection in inclination 
and eccentricity spaces. 



fl2008h . These 



Figure 10 allows direct comparison to the results of iRoig et al. 
authors calculate size- frequency distributions for the possible number of Vestoids 
that have migrated across the 3:1 resonance into dynamically stable orbits (see 
their Fig. 6). In their most optimistic simulations they find that on the order 
of 100 Vestoids with a size of 2 km or greater could cross this resonance over 
the maximum lifetime of the Vesta family (3.5 Gyr). Our lower and upper 
limits suggest that between 140 and 2800 basaltic asteroids with a size bigger 
than 2 km (if < 15.0) exi st beyond 2 . 5 AU . The lower end of this estimate 



is commensurate with the 



Roig et al 



(I2OO8I ) result. However, if the basaltic 



asteroids in the outer belt are exclusively of Vestian origin then they would 
be distributed at semi- major axes near the 3:1 resonance. Our calculated dis- 
tribution is in contrast to this expectation: it is roughly uniform across all 
a > 2.5 AU and has a single peak at 2.7 AU. This peak might be related 
to basaltic ast eroids generated from a partially differe ntiated Eunomia parent 



body (see §6.1, 



Reed at el 



1997 



Nathues et al 



2OO5I ) . Ignoring this peak, the 



flat "background" population beyond 2.5 AU is suggestive of multiple sources 
of basaltic asteroids in addition to those few added by the diffusion of Vestoids 
across the 3:1 resonance. The very low number statistics in this region make 
it difficult to make more detailed inferences. Further observations are needed 
to provide more robust constraints on the distribution of these objects beyond 
the 3:1 resonance. 
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6 Discussion and Summary 



We have identified 50 candidate basaltic asteroids from tiie SDSS MOC based 
on pliotometric and dynamical constraints. Based on observations of these 
basaltic candidates we have constrained the amount of basaltic material that 
is present outside of a region occupied by the dynamical Vestoids (defined by 
a-e-i orbital elements, §2). We have determined the unbiased distribution of 
this basaltic material in absolute magnitude and semi-major axis. 



6.1 Basaltic Asteroids with a > 2.5 AU 



The prospect of multiple basaltic asteroids with semi-major axes greater than 
2.5 AU is particularly compelling due to the unlikelihood of numerous such 
objects originating from the surface of Vesta and then cro ssing the 3:1 mean 



motion resonance (iNesvorny et al 



2008 



Roig et al 



20081 ) 



Our selection process identified 12 basaltic candidates between the 3:1 and 
5:2 mean motion resonances at 2.5 and 2.8 AU respectively. Six of these are 



clustered around the Eun omia family 



families in the Main Belt (jZappala et al. 



ig. 3) which is one of the largest 



19951 ) and may have derived from 



a partially or fully differentiated parent body. This clair n is suppor t ed by 
a compositional gra dient across the surface of Eunomia (IReed at el.l . 119971 : 



Nathues et al. 



20051 ) a nd a spread of surfac e compositions amongst members 



of the Eunomia family (ILazzaro et al 



a differentiated mineralogy ( 



Gaffey et al 



19991) . b oth of which are consistent with 



2OO2I ) . Presumably the basaltic crust 



of the Eunomia parent body was removed in the large family forming collision 
and/or by subsequent collisions. Asteroid 21238 (1995 WV7) is a spectroscop- 
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i cally confirmed V- 



t ype that may represent some of this los t basa ltic crust 



( iHammergren et al. 



2006 



Binzel et al. 



2007 



Carruba et al. 



2003). One of 



the six candidates in the vicinity of Eunomia, asteroid 40521, is a spectroscop 
ically confirmed V-type, however there is some uncertainty abo ut whether this 
object originated from the Vesta or Eunomia parent bodies fiCarruba et al. 



20071 : 



Roig et al. 



20081 ). Detailed numerical simulations will be necessary to 
investigate whether it is dynamically feasible for the remaining five candidates 
to have originated from Eunomia. 

It is very unlikely that the selection of 6 candidates in the Eunomia family 
represents a statistical fluke related to the high density of objects in this region. 
Approximately 11% of MOC objects beyond the 3:1 resonance are part of the 
Eunomia family. Thus, assuming a binomial probability distribution for which 
there is an 11% chance of a single candidate being selected at random from the 
Eunomia family, there is only a 0.7% probability that 6 of the 16 candidates 
from beyond the 3:1 resonance would belong to the Eunomia family. 

We have identified only 4 basaltic candidates with semi-major axes greater 
than 2.8 AU (Fig. 3, Table 2). One of these, asteroid 10537 (1991 RY16), may 
be classified as a V-type because of its deep 1 fim absorption feature which 
is typical of this class and indicative of a pyroxene- rich and/or olivine-rich 
mineralogy. However it also has an unusually deep absorption feature cen- 
tered at 0.63 fim that is uncommon to V-types (Fig. 5). These visible wave- 
length data have been confirmed by other observations [Duffard & Roig (2008), 
submitted to Planetary and Space Science]. Follow-up infrared observations 



provide additiona . 



( jMoskovitz et al. 



insight on the surface composition of this unusual object 



20081 ). Another outer belt candidate is unnumbered asteroid 



2003 SG55. It is possible that this object originated from a differentiated Eos 
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parent body (iMothe-Diniz and Carvano . 



20051 ). The pre sence of multiple non 



linear seen 



1990 



ar resonances that pass through the region (iMilani and Knezevic 



1992) may have provided a 



SG55 even though HCM analysis (IZappala et al. 



pathway for the or 



Dital migration of 2003 



1990) shows that 2003 SG55 



is outside of the Eos dynamical family. Future study of these outer belt basaltic 
candidates will have important implications for understanding the dynamical 
and thermal histories of these specific families. 



6.2 Basaltic Asteroids with a < 2.5 AU 



Defining the dynamical boundaries of the Vestoi ds is difficult due to effect s 



2005 



2003). 



such as Yarkovsky drift and resonance scattering (iCarruba et al. , 
Nevertheless we have attempted to identify non-Vestoid basaltic candidates in- 
side of 2.5 AU because neglecting to do so could exclude basaltic asteroids of 
non-Vestian origin. The relatively uniform distribution of basaltic candidates 
in this region suggests two possibilities: (1) there exists non-Vestoid basaltic 
asteroids in the inner belt that represent a unique differentiated parent body, 
or (2) Vestian fragments actually populate the full dynamical extent of the 



inner Main B elt. The first possibi 



analyses (e.g. 



Lawrence and Luceyl (120071 ): 



i ty is testable by performin 



Gaffey et al. 



g mineralogical 



fll993ah ) on NIR spec- 



tra of basaltic candidates with large values of Av, for example 38070 (1999 
GG2) with Av = 2.1 km/s and 33881 (2000 JK66) with Av = 3.3 km/s. If 
these inner belt basaltic candidates do represent a unique differentiated par- 
ent body then they may be compositionally distinguishable from Vesta and 



the Vestoids (e.g. Magnya, 



Hardersen et al. 



(12004)'). The second poss ibility 



can be investigated with numerical simulations. iNesvorny et al. 



(120081) show 



33 



that the dynamical diffusion of Vestian material does in fact populate the 
semi-major axis extent of the inner Main Belt, but does not extend below 
orbital inclinations of 4°. In the future, the most effective test for inner belt, 
non-Vestoid basaltic asteroids will be a combination of dynamical studies and 
mineralogical analyses. 



6.3 Total Basaltic Asteroid Inventory 



We find that the calculated size-frequency distribution of basaltic asteroids 
outside of the dynamical Vestoid population is heavily weighted towards small 
bodies. We have also shown that the mass of basaltic material outside of the 
dynamical Vestoid population is M^ur < 5.9 x 10^^ kg. This value is comparable 
to the mass of basaltic material inside of the dynamical Vestoid population. 
Due to uncertainties in the definition of dynamical families we do not attempt 
to make a statement about the fraction of M^ur that is truly non-Vestoid 
in origin. We have attempted to place upper limits on Msur and My and a 
conservative estimate on M^x- Even so we note that Ms^^-I-My < M^x by over 
an order of magnitude. 



The size-frequency distribution {N{H)) and smaller total mass relative to Mex, 
both inclusive and exclusive to the Vestoid dynamical family, are consistent 
with most (but not all!) of the basaltic asteroids derivin g from the Ves t a par- 



19971 ) 



ent body. Vesta has several large craters on its surface (iThomas et al. 
whose coUisional formation would have produced fragments of varying size. 
Due to ejection velocities, the Yarkovsky effect, resonance interactions and 
close encounters, these fragments d i ffused away from the Vesta parent body 



(jBottke et al. 



2001 



Carruba et al. 



2005 



Nesvorny et al. 



20081 ). The small- 
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est fragments experienced the fastest dynamical evolution while larger objects 
would take longer to separate from Vesta due to t he decreased s i gnific ance of 



2000). Over 



the Yarkovsky force for bodies larger than 10 km (iBottke et al.l . 
time, small Vestoids migrated to the edge of the family at the z/g and 3:1 res- 
onances. Although many of the fragments that reached these resonances were 



scattered into planet crossing orbit s or into th e Sun^ some crossed over int o 



stable orbits (iGladman et al. 



1991 



Roig et al. 



2008 



Nesvorny et al. 



200811. 



Thus, the size-frequency distribution of Vestoids shows a depletion of small 
bodies {H > 16) relative to the SDSS MOC distribution. This is apparent in 
Fig. 9 where the slope of the Vestoid distribution (solid curve) for if > 16 
is slightly steeper than that of the MOC (dotted curve). This also holds true 
when compared to just the subset of MOC objects that are confined to the 
semi-major axis range of the Vestoids (2.2 < a < 2.5 AU). 

The overabundance of small bodies (sizes < 2 km) amongst the basaltic aster- 
oids outside of the dynamical Vestoid population (Fig. 7) may represent the 
hiding place of these "missing Vestoids". If this is true then it suggests that 
the shallow slope of absolute ma gnitude distributions at large values of H for 



HCM-defined collisional families (iMorbidelli et al. 



20031 ) may be attributed to 



the preferential loss of small bodies out of the families by dynamical processes 
such as the Yarkovsky effect and resonance scattering. 

The semi-major axis distribution indicates that the full inventory of basaltic 
material in the Main Belt is not consistent with an origin exclusively cou- 
pled to the formation of the large Vestia n crater. With an age of 1-3.5 Gyr 



( jBottke et al. 



2005a 



Marzari et al. 



19961 ). the primary Vesta crater was formed 



after the completion of planet formation. Thus, the asteroidal fragments gener- 
ated in this impact evolved in a dynamical environment similar to the present 
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day architecture of secular and mean motion resonances. In such an environ- 
ment it is relatively easy for asteroids to migrate in semi-major axis because of 
the Yarkovsky effect, however it is more difficult to produce significant changes 
in inclination and eccentricity. This motivates the search for basaltic asteroids 
at very different i and e values than Vesta. 



Our results suggest that basaltic asteroids with si zes in exces s of 2 



im do exist 



Roig et al. 



(120081 ) show that 



at low inclinations and beyond 2.5 AU (Fig. 10). 
it is not possible to dynamically transport many hundreds of Vestoids with di- 
am eters greater than 2 krn beyond the 3:1 resonance. This is further confirmed 



by 



Nesvorny et al. 



(120081 ) whose dynamical simulations do not produce any 
3 km size bodies beyond this resonance or any Vestoids with low inclination 
orbits a < 4°). 



We cite three explanations for the presence of basaltic material in regions of 
phase space that are dynamically inaccessible to the Vestoids. One possibility 
is that these basaltic asteroids originated from the surface of Vesta, not from its 
largest crater but sometime before the Late Heavy Bombardment (LHB) more 
that 3.8 Gyr ago. The higher density of objects in the pre-LHB asteroid belt 



(iPetit et al. 



2OO2I ) would have made for a dynamical environment in which 
coUisional and scattering events were greatly enhanced over present day rates. 
In addition, resonance sweeping during the LHB could have dispersed these 



asteroids' inclinations (iNesvorny et al 



20081) 



A second possibihty is that these objects represent sc attered planetesima. 



20061 ). 



from the terrestrial planet region of the Solar System (iBottke et al. 
In the inner Solar System solid-body accretion times would have been faster 
than in the Main Belt. These bodies would have accreted a higher quantity 
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of live radioactive nuclides such as ^^Al and been more likely to differentiate. 
The amount of basaltic material that was injected into the Main Belt from 
the inner Solar System is currently unconstrained. 

A final explanation for the presence of basaltic asteroids at low inclinations 
and beyond 2.5 AU is that they derived from other differentiated parent bodies 
that formed in the Main Belt. Although we find no asteroidal evidence for a 
large number of undiscovered basaltic asteroids, we cannot rule out minor 
contributions to the basaltic inventory from non-Vestoid differentiated bodies 



whose basalt has been heavily dep^ 
(i.e. a "battered to bits" scenario 



eted due to col 



Burbine et al. 



isions and dynamical effects 



19961 ). It is possible that the 



observed distribution of basaltic asteroids in the Main Belt is derived from 
many different source populations including the Vestoids, inner Solar System 
planetesimals and a heavily depleted population of non-Vestoid differentiated 
bodies. 



6.4 Comparison to Other Works 



We comp are our calculated sem i-ma jor axis distr i bution to the simulation 



results of 



Nesvorny et al. 



(120081 ) and 



Bottke et al. 



(120061 ) (Fig. 11). This al- 



lows us to investigate the feasibility of two different scenarios for delivering 
basaltic asteroids to orbits outside of the dynamical Vestoid population. Fig- 
ure 11 shows the normalized distribution of non-Vestoid basaltic asteroids as 
a function of semi-major axis. These distributions have been normalized by 
their maximum values (all at 2.2 AU) because limitations inherent to each 
work precludes the comparison of an absolute number of basaltic asteroids. 
We confine these distributions to the inner Main Belt (a < 2.5 AU) because of 
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l ow nu mber statistics for larger values of a in our analysis and the lBottke et al. 



the 



(120061) simulations and th e lack of any Vestoids across the 3:1 resonance in 



Nesvorny et al. 



(120081 ) simulations. The distributions from the numerical 



simulations exclude test particles that populate the dynamical Vestoid region 
of orbital element space as defined in §2. This produces distributions that are 
non- Vestoid based upon our specific dynamical criteria. 



The iNesvorny et al.l (120081 ) study was designed to understand the dynamical 
evolution of Vesta's ejected fragments. They modeled the diffusion of Vestoids 
over an estimated family lifetime of 2 Gyr. If our distribution of basaltic 
asteroids was exclusively deri ved from dynam i cally a berrant Vestoids then we 
would expect it to match the 



Nesvorny et al. 



mm result. 



The lBottke et al.l (120061 ) study models the injection of planetesimal-size bodies 
from the inner Solar System into the Main Belt during the epoch of planet 
formation. The specific source region of these planetesimals determines which 
are differentiated and contain basalt. Due to uncertainties in the time-scales 
of accretion and the initial abundance of radioactive isotopes, it is not clear 
where the boundary between differentiated and u ndifferentiated plane tesimals 



Bottke et al. 



(120061 ) results: 



existed. Thus we plot two distributions from the 
one corresponding to all of the planetesimals that are scattered into the main 
belt from initial semi-major axes, oq, less than 1.5 AU and those with ao < 2.0 



AU. If our distribution was deriv ed primarily 



we would expect it to match the iBottke et al. 



rom t 



lese scattered interlopers 



mm result. 



We calculate the chi-squared (x^) probability tha t these simulated dist ribu 



tions are derived from our observational result. The 



tribution and the .Bottke et al. 



Nesvorny et al. 



mm dis- 



(120061) distribution for a planetesimal source 
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region of < 1.5 AU produce probabilities equal to 99.6% and 98.5% 
respectively. Within the limits of our calculations these distrib utions fit the 



data e qually well. The worst fit (x^ probability = 73%) is the 



Bottke et al. 



( 120061 ) distribution for a planet esimal source 



ference in probabilities from the 



Bottke et al 



region of ao < 2.0 AU. The dif- 



(120061 ) simulations suggests that 



the contribution of basaltic asteroids from inner Solar System planetesimals 
is weighted towards objects originally separated in semi-major axis from the 
Main Belt rather than objects that formed at semi-major axes around 2.0 AU. 
This is not unexpected as faster accretion times at smaller heliocentric dis- 
tances would have produced a greater number of differentiated planetesimals. 



This work is the first survey specifically designed to constrain the global dis- 
tribution of b asaltic material in t he Main Bel t. Future all-sky sur veys such as 



PanSTARRS (iKaiser et al 



200J) will greatly 



120021 ) and LSST (IClaver et all 
increase the number of known asteroids and will allow a more complete de- 
scription of the basaltic asteroids. The results presented here suggest that the 
total amount of basaltic material in the Main Belt is very much less than 
expected from the ~ 100 differentiated parent bodies traced by the iron me- 
teorites and basaltic achondrites, a result which agrees with the suggestion 
that basalti c fragments have bee n comminuted to sizes below observational 



limits (e.g. IBurbine et al 



19961 ) or that some unknown space weathering 



processes has obscured t h e spe ctral signature of basalt on asteroid surfaces 



(jWetherill and Chapman 



19881 ). However, the inferred presence of a small 
amount of basaltic material throughout the Main Belt suggests that some 
vestige of crustal material from large differentiated bodies does exist and may 



represent the upper size limit of this population. 
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Color Range 



u- g 


[0.170,0.350] 


g-r 


[0.201,0.281] 


r — i 


[0.061,0.121] 


i — z 


[-0.389,-0.249] 


PCI 


[0.360,0.520] 


PC2 


[-0.258,-0.078] 



Table 1 

Vestoid Colors - 50*'* Percentile 
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Number Designation a (AU) e sin(i) Av (km/s) RG06? H06? B07? 

156914 2003 FB31 2.37 0.10 0.14 0.3 

28517 2000 DD7 2.29 0.09 0.14 0.4 / / / 

124343 2001 QF117 2.22 0.11 0.12 0.7 / 

112257 2002 LC13 2.22 0.10 0.11 0.7 / 

54668 2000 W085 2.44 0.11 0.10 0.7 / 

10544 1992 DA9 2.23 0.11 0.13 0.8 / / 

7558 Yurlov 2.29 0.11 0.09 0.9 / 

74894 1999 TB119 2.22 0.12 0.14 0.9 / 

67416 2000 QW64 2.17 0.08 0.10 1.0 / 

65507 4151 P-L 2.22 0.12 0.11 1.0 / / / 

60669 2000 GE4 2.21 0.13 0.12 1.0 / / / 

24941 1997 JM14 2.48 0.12 0.09 1.0 / 

55804 1994 PD13 2.29 0.13 0.15 1.1 

2823 van der Laan 2.41 0.06 0.08 1.1 / / / 

56570 2000 JA21 2.38 0.10 0.07 1.2 / / / 

98278 2000 SE212 2.16 0.06 0.10 1.3 / / 
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55998 1998 SQ135 2.19 0.10 0.08 1.3 / / 

46690 1997 AN23 2.24 0.10 0.07 1.3 / / 

2002 RS123 2.37 0.16 0.14 1.4 / / 

6081 Cloutis 2.49 0.16 0.12 1.4 / / 

110544 2001 TX96 2.68 0.06 0.17 1.6 

10537 1991 RY16 2.85 0.10 0.11 1.6 / / 

122656 2000 RV93 2.75 0.13 0.09 1.7 / 

44114 1998 HN21 2.18 0.02 0.10 1.7 / 

8805 1981 UMll 2.34 0.12 0.05 1.7 / / 

2002 TF149 2.77 0.02 0.10 1.9 / 

109387 2001 QL168 2.59 0.13 0.20 2.0 / 

86493 2000 DD17 2.14 0.09 0.05 2.0 

158606 2002 XS80 2.26 0.17 0.07 2.1 

38070 1999 GG2 2.14 0.13 0.06 2.1 

2001 QA244 2.73 0.18 0.17 2.2 

93322 2000 SA221 2.74 0.10 0.21 2.2 / 

84036 2002 PF50 2.33 0.15 0.04 2.2 / / 
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55092 2001 Q0123 

69634 1998 FH68 

41910 2000 WS141 

40521 1999 RL95 

152169 2003 AV27 

107008 2000 YM112 

2003 SG55 

50802 2000 FH27 

140706 2001 UV79 

55550 2001 XW70 

172731 2004 BY120 

129632 1998 HV36 

84021 2002 PC41 

33881 2000 JK66 

86627 2000 EJ126 

2001 GQ5 

111515 2001 YC91 



2.15 0.13 0.05 
2.38 0.16 0.05 
2.26 0.20 0.13 
2.53 0.05 0.22 
2.25 0.10 0.02 
2.35 0.18 0.05 
3.01 0.13 0.20 
2.91 0.09 0.03 
2.61 0.03 0.02 
2.55 0.14 0.23 
2.74 0.14 0.23 
2.53 0.24 0.12 
2.68 0.09 0.25 
2.21 0.24 0.18 
2.42 0.24 0.04 
2.20 0.26 0.15 

3.16 0.16 0.28 



2.2 / 

2.3 / 
2.3 / / 

2.3 / / 
2.4 

2.4 / 

2.5 / 

2.5 / 

2.6 / 
2.6 / 
2.8 

3.0 

3.0 / 

3.3 / / 
3.5 
3.6 

3.9 / 
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Table 2: Non-Vestoid Basaltic Candidates 
The table columns are: asteroid number and designa- 
tion, semi-major axis, eccentricity, sine of inclination, 
ejectio n velocity (Af ) requ ired to have originated from 



Vesta (IZappala et al. 



19961 ) and the last three columns 



indicate whether the asteroid has been identified as a 



V-type candidate b v iRoig and Gil-HuttonI ( 



Hammergren et al 



(12006 



20061. RG06) 



(12007 . 



H06) and iBinzel et al 
B07). The candidates are sorted by their Av values. Spec- 
troscopic data have been obtained for candidates listed in 
bold- face (Table 3, Fig. 4). The orbital elements of Vesta 
are: a = 2.36 AU, e = 0.09 and sm{i) = 0.12. 
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Object 



UT Date Inst. Mag. texp (s) Analog Type 



7558 (Yurlov) 



Oct. 1, 2006 ESI 16.1 1800 Ave 



V 



10537 (1991 RY16) Oct. 1, 2006 ESI 17.9 2700 Avg. V? 



Jan. 18, 2007 SNIFS 18.3 3600 SA93-101 V? 



24941 (1997 JM14) Oct. 1, 2006 ESI 19.0 900 Avg. V 



28517 (2000 DD7) Mar. 10, 2007 ESI 16.5 900 Avg. V 



33881 (2000 JK66) Nov. 23, 2007 SpeX 16.6 2880 HD9562 V 



38070 (1999 GG2) Aug. 23, 2006 ESI 17.8 2700 Avg. V 



Oct. 05, 2006 SpeX 16.9 1920 HD377 V 



46690 (1997 AN23) Mar. 10, 2007 ESI 18.7 2700 Avg. 



56570 (2000 JA21) Jan. 4, 2006 ESI 19.0 2700 Avg. 



V 



60669 (2000 GE4) Jun. 18, 2007 SNIFS 17.8 4500 BS5534 V 
2823 (van der Laan)^ Nov. 22, 2005 SpeX 16.4 480 SA93-101 V 



40521 (1999 RL95)2 Apr. 30, 2006 CMOS 18.1 3000 SA107-871 V 
Table 3 

Observation Summary 

The columns in this table are: object number and designation, UT date of obser- 
vation, the instrument used, magnitude of the target from the JPL HORIZONS 
system, net exposure time in seconds, the solar analog that was used for calibration 
in the reduction process (see §4, Avg. indicates that the average of all solar analogs 
for a given night were used) and the taxono mic designation t hat w e assigned to the 



spectra based upon the criteria outlined by 



Bus and Binzell (12002) 



^SMASS Database ( |http: / /smass.mit .ecjiy^minus.htrnll ) , IRTF 
Roig et al.1 (120081), Gemini South 



Object 



H Type Reference 



63 (Ausonia) 


7.25 


S 


(1) 


556 (Phyllis) 


9.23 


S 


(1) 


1145 (Robelmonte) 


11.05 


TDS 


(2) 


2098 (Zyskm) 


11.92 


TX 


(2) 


2086 (Newell) 


11.99 


Xc 


(1) 


2346 (Lilio) 


12.04 


C 


(1) 


3376 (Armandhammer) 


12.33 


Sq 


(1) 


3865 (1988 AY4) 


12.42 


Xc 


(1) 


4510 (Shawna) 


12.69 


s 


(2) 


2024 (McLaughlin) 


12.73 


s 


/ <r\\ 

(2) 


1781 (van Biesbroeck) 


12.75 


xs 


^ <r\\ 

(2) 


5111 (Jacliff) 


12.78 


R 


(1) 


4845 (Tsubetsu) 


12.81 


X 


(1) 


1697 (Koskenniemi) 


12.93 


TX 


(2) 


5600 (1991 UY) 


12.94 


S 


(3) 



Table 4 

Int erlopers in Vesto i d Dy namical Space 



(1) 
(2) 
like: 

(3) 



3us and Binzel 



Xuetal 



2002) 



19951 ). uncertain classifications are listed in order of taxonomic 



ihood 



Lazzaro et al 



(|2004l ) 
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Fig. 1. Optical (ESI on Keck II) and near-infrared (SpcX on NASA's Infrared 
Telescope Facility) spectra of basaltic asteroids 38070 (1999 GG2) and 33881 (2000 
JK66) and a laboratory spectrum of polymict eucrite Y74450 (acquired at the NASA 
RELAB facility, Brown University by Takahiro Hiroi) . These data clearly show the 
broad 1 and 2 fj,m absorption features that are characteristic of basaltic material. 
The band centers of the SDSS ugriz filters are indicated on the spectrum of 38070. 
The spectra have been offset by increments of ±0.5 units relative to the 1999 GG2 
spectrum. 
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Fig. 2. i — z color histogram of the 1051 dynamical Vestoids that are included in 
the SDSS MOC. A gaussian fit to the distribution is shown by the dashed curve. 
The range of colors that correspond to the 50*'* percentile of dynamical Vestoids 
(defined symmetrically relative to the center of the gaussian) is indicated. 
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Semi — major Axis (AU) 



Fig. 3. Osculating orbital elements of Vestoids (grey circles), MOC objects (small 
dots) and basaltic candidates (x). The Eunomia and Eos families are indicated. 
The apparent overlap of some candidates with the Vestoids is due to projection 
along the unplotted eccentricity axis. Note that there are at least 6 candidates in 
the vicinity of the Eunomia family. 
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Fig. 4. Optical spectra of observed basaltic candidates. Each spectrum has been 
normalized at 0.55 /Ltm. Due to instrumental problems, only the red channel data 
from SNIPS is shown for asteroid 60669 (2000 GE4). A spectrum of Vestoid 6093 
Makoto is included for reference (obtained with SNIPS on 26 October 2006). The 
shallow depth of the 1 iim band in the spectrum of 1997 AN23 suggests that it is 
an S-type asteroid, whereas all of the others are V- types. 
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0.4 0.5 0.6 0.7 0.8 0.9 

Wavelength (/Ltm) 

Fig. 5. ESI and SNIFS spectra of basaltic candidate asteroid 10537 (1991 RY16). 
The gap in the SNIFS data corresponds to the transition region between the red 
and blue channels. The SNIFS data have been offset for clarity. 
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102 I 

13 14 15 16 17 

H 

Fig. 6. Differential number distributions as a function of absolute magnitude for the 
SDSS MOC and the ASTORB database. The x-axis has been limited to the range 
of absolute magnitudes represented by our basaltic candidates. The extrapolation 
of Eq. 4 is indicated by the dashed line. The ratio of the SDSS MOC distribution 
to Naii{H) gives the completeness and efficiency of the SDSS as a function of H. 
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13 14 15 16 17 

Absolute Mognitude 

Fig. 7. Unbiased differential number distribution as a function of absolute magni- 
tude for basaltic asteroids in the Main Belt, excluding those within the dynamical 
Vestoid population. 90% Poissonian confidence levels are plotted as upper and lower 
limits for each point. Certain H-h'ms are empty because they contained no candi- 
dates. Upper limits are indicated by arrows for these bins. This distribution is well 
represented by a single function of the form k10°'^ (dotted line). 
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102 I 

13 14 15 16 17 

Absolute Magnitude 

Fig. 8. Unbiased cumulative number distribution as a function of absolute magni- 
tude for basaltic asteroids in the Main Belt. The value for each bin represents an 
upper limit to the number of asteroids within the respective absolute magnitude 
range, thus producing a liberal approximation to the cumulative number. 
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Fig. 9. Number distribution of dynamically defined Vestoids as a function of abso- 
lute magnitude (solid). The fit to the distribution for 13.4 < H < 14.8 is indicated 
by the dashed line. The dotted curve is the SDSS MOC distribution. Note the in- 
creased deficiency of small bodies {H > 16) in the Vestoid family relative to the 
MOC. 
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Fig. 10. Number of basaltic asteroids as a function of semi-major axis and abso- 
lute magnitude. The integrated number of objects has been scaled to equal those 
represented by Fig. 8. Three different absolute magnitude regimes are shown. The 
upper and lower limits are calculated from 90% Poisson confidence levels. For bins 
without any candidates N{a; H) becomes zero. Upper limits are indicated with ar- 
rows for these bins. Assuming a typical Vestian albedo of 0.4, the three absolute 
magnitude limits correspond to si zes of 0.76 km, 1 . 3 km and 2.0 km (from top to 



bottom). Dynamical simulations (jRoig et al 



20081 : 



Nesvornv et al. 



2008|) predict 



that Vestoids smaller than 3 km can cross the 3:1 resonance at 2.5 AU. 
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Fig. 11. Comparison of semi-major axis distributions. Our calculated distribution 
of basaltic asteroids is shown as the solid line with 90% Poissonian upper and lower 
limits. The results of two dynamical simulations are over-plotted. The dotted line 
shows the simulated number of Vestian fragme nts that have evolved outside of the 



2003). The dashed 



HCM-defined dynamical limits of this family (jNesvornv et al. 
lines plot the distribution of inner Solar System planetesimals scattered into the 
Main Belt from two different s ource regions , a < 1.5 AU (black dashed line) and 



a < 2.0 AU (grey dashed hue) (jBottke et al 



2009 1 
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